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(Dated: April 30, 2007) 

We have measured the Ab lifetime using the exclusive decay Ab — ► J /if) A, based on 1.2 fb _1 of 
data collected with the DO detector during 2002-2006. From 171 reconstructed A b decays, where 
the J/tp and A are identified via the decays J/tp — > and A — > pix, we measured the Ab lifetime 

to be t(A 6 ) = 1.218i° "g(stat) ± 0.042(syst) ps. We also measured the B° lifetime in the decay 
B° -f J/iP(n + Li-)K%(-K + ir-) to be t(B°) = 1.50llo;° 7 8 (stat) ± 0.050(syst) ps, yielding a lifetime 
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ratio of r(A 6 )/r(B°) = 0.811±^ (stat) ± 0.034(syst). 
PACS numbers: 14.20.Mr, 14.40.Nd, 13.30.Eg, 13.25.Hw 



Lifetime measurements of b hadrons provide important 
information on the interactions between heavy and light 
quarks. At leading order in Heavy Quark Effective The- 
ory (HQET) [lj], light quarks are considered spectators 
and all b hadrons have the same lifetime. Differences 
arise at higher orders when corrections from interactions 
are taken into account. For HQET calculations of order 
l/m^, where mj is the mass of the b quark, the agree- 
ment between the predicted lifetimes and the experimen- 
tal results is excellent for B mesons |2j. However, in the 
b baryon sector, the world average of measurements of 
r(A 6 )/r(B°) = 0.844 ± 0.043 3 is smaller than the pre- 
diction of the ratio at this order. Recently there have 
been significant improvements in theoretical calculations 
ofr(A 6 )/r(S ). Next-to-leading order effects in QCD 0], 
corrections at (9(l/m£) in HQET and lattice QCD 
studies [a], have led to an improved theoretical predic- 
tion, r(A 6 )/r(S ) = 0.88 ±0.05 7]. This value agrees 
with previous experiments to within the current theoret- 
ical and experimental uncertainties. However, a recent 
precise measurement Q reports a value of the Af, life- 
time consistent with B meson lifetimes, and the ratio 
r(A;,)/r(i3 ) consistent with unity. Additional precise 
measurements of the Af, lifetime and t(A{,)/t(B ) ratio 
may help settle this question. 

In this Letter, we report measurements of the Af, life- 
time using the exclusive decay At, — > J/tp + A, and its 
ratio to the B° lifetime using the B° — ► J/ip + Kg de- 
cay channel. This B° decay channel is chosen because of 
its similar topology to the Af, decay. The J/ip is recon- 
structed from the n + fJ-~ decay mode, the A from ptt~ , 
and the Kg from tt + it~ . Throughout this Letter, the ap- 
pearance of a specific charge state also implies its charge 
conjugate. The data used in this analysis were collected 
during 2002-2006 with the DO detector in Run II of the 
Tevatron Collider at a center-of-mass energy of 1.96 TeV, 
and correspond to an integrated luminosity of 1.2 fb _1 . 

The DO detector is described in detail elsewhere Q. 
The detector components most relevant to this analysis 
are the central tracking and the muon systems. The for- 
mer consists of a silicon microstrip tracker (SMT) and a 
central scintillating fiber tracker (CFT) surrounded by a 
2 T superconducting solenoidal magnet. The SMT has a 
design optimized for tracking and vertexing for pseudo- 
rapidity of |ry| < 3 [icj |. For charged particles, the reso- 
lution on the distance of closest approach as provided by 
the tracking system is approximately 50 /jm for tracks 
with pt ~ 1 GeV/c, where px is the component of the 
momentum perpendicular to the beam axis. It improves 
asymptotically to 15 fim for tracks with pt > 10 GeV/c. 
Preshower detectors and electromagnetic and hadronic 



calorimeters surround the tracker. The muon system 
is located beyond the calorimeter, and consists of mul- 
tilayer drift chambers and scintillation counters inside 
1.8 T toroidal magnets, and two similar layers outside 
the toroids. Muon identification for \r/\ < 1 relies on 10 
cm wide drift tubes, while 1 cm mini-drift tubes are used 
for 1< \rj\ < 2. 

The primary vertex of the pp interaction is determined 
for each event using the average position of the beam- 
collision in the plane perpendicular to the beam as a 
constraint. The precision of the primary vertex recon- 
struction is on average 20 /im in the plane perpendicular 
to the beam and about 40 fj,m along the direction of the 
beam. 

We base our data selection on reconstructed charged 
tracks and identified muons. Although we do not require 
any specific trigger, most of the selected events satisfy 
dimuon or muon triggers. To avoid a trigger bias in 
the lifetime measurement, we reject events that depend 
on impact parameter based triggers. We start the Af, 
and B° reconstruction by searching for events with J/ %p 
mesons. We then search in these events for A and Kg 
particles. To reconstruct J/ip — > candidates, we 

select events with at least two muons of opposite charge 
reconstructed in the tracker and the muon system. The 
track of each muon candidate must either match hits in 
the muon system, or have calorimeter energies consistent 
with a minimum-ionizing particle along the direction of 
hits extrapolated from the tracking layers. For at least 
one of the muons, we require hits in all three layers of 
the muon detector. Both muons are required to have 
Pt >2.5 GeV/c if they are in the region \rj\ <1. The 
muon tracks are constrained to originate from a com- 
mon vertex with a x 2 probability greater than 1%, and 
each J /%p candidate is required to have a mass in the 
range 2.80-3.35 GeV/c 2 . The A — > pir~ decays are recon- 
structed from two tracks of opposite charge constrained 
to a common vertex with a x 2 probability greater than 
1%. Each A candidate is required to have a mass in the 
range 1.100-1.128 GeV/c 2 . The proton mass is assigned 
to the track of higher px, as observed in Monte Carlo 
studies. To suppress contamination from cascade decays 
of more massive baryons such as E° — > A7 or S° — > Air , 
we require the cosine of the angle between the px vec- 
tor of the A and the vector in the perpendicular plane 
from the J/tp vertex to the A decay vertex to be larger 
than 0.9999. For A's that decay from Af, the cosine of 
this angle is very close to one. The Kg — > it + tt~ selec- 
tion follows the same criteria, except that for the Kg, 
the mass window is 0.460-0.525 GeV/c 2 , and pion mass 
assignments are used. 
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We reconstruct the A& and B° by performing a con- 
strained fit to a common vertex for either the A or K s 
and the two muon tracks, with the latter constrained 
to the J/i[> mass of 3.097 GeV/c 2 Q. Because of their 
long decay lengths, a significant fraction of A and K s 
particles will decay outside the SMT. There is therefore 
no requirement of SMT hits on the tracks from A and 
K s decays. To reconstruct the Af, (B°), we first find 
the A (Kg) decay vertex, and then extrapolate the mo- 
mentum vector of the ensuing particle and form a vertex 
together with the two muon tracks belonging to the J/ip- 
If more than one candidate is found in the event, the can- 
didate with the best \ 2 probability is selected as the A;, 
(B°) candidate. The mass is required to be within the 
range 5.1-6.1 GeV/c 2 for Af, candidates and within 4.9- 
5.7 GeV/c 2 for B° candidates. For the choice of the final 
selection criteria, we optimize S/ \/S + B, where S and 
B are the number of signal (A&) and background candi- 
dates, respectively, by using Monte Carlo estimates for^S* 
and data for B. For the Monte Carlo, we use PYTHIA [11 1 
and EVTGEN [l2j to produce and decay particles, respec- 
tively, and GEANT3 [l3[ to simulate detector effects. As 
a result of this optimization, the px of the A (K s ) is re- 
quired to be greater than 2.4 (1.8) GeV/c, and the total 
momentum for both Af, and B° is required to be greater 
than 5 GeV/c. Finally, any candidate which has been 
identified as a Af, is removed from the B° sample. 

We determine the decay time of a Af, or B° by mea- 
suring the distance traveled by the b hadron candidate in 
a plane transverse to the beam direction, and then ap- 
plying a correction for the Lorentz boost. We define the 
transverse decay length as L xy = L xy ■ pt /pt where L xy 
is the vector that points from the primary vertex to the 
b hadron decay vertex and pt is the transverse momen- 
tum vector of the b hadron. The cvent-by-event value of 
the proper transverse decay length, A, for the b hadron 
candidate is given by: 



A 



L X y _ j CM B 



(1) 



where (/?7)y an d Mb are the transverse boost and the 
mass of the b hadron. In our measurement, the value of 
M B in Eq. [1] is set to the Particle Data Group (PDG) 
mass value of A& or B° Q. We require the uncertainty 
on A to be less than 500 /xm. 

We perform a simultaneous unbinned maximum like- 
lihood fit to the mass and proper decay length distribu- 
tions. The likelihood function L is defined by: 



C = 



n 



N 



■ exp (~n s - n b ) 



b -rJ 

hn b bkg 



(2) 



where n s and rif, are the expected number of signal and 
background events in the sample, respectively. N is the 
total number of events. T 3 sig (-T^kg) i s the product of 



three probability density functions that model the mass, 
proper decay length, and uncertainty on proper decay 
length distributions for the signal (background). We di- 
vide the background into two categories, prompt and 
non-prompt. The prompt background is primarily due 
to direct production of J ftp's which are then randomly 
combined with a A or Kg candidate in the event. The 
non-prompt background is mainly produced by the com- 
bination of J/ip mesons from b hadron decays with A or 
Kg candidates present in the event. 

For the signal, the mass distribution is modeled by a 
Gaussian function, and the A distribution is parametrized 
by an exponential decay convoluted with the resolution 
function: 



G(X j ,a j 



1 



2-Ksa. 



■ exp 



-\ 2 



2(sa 3 Y 



(3) 



where \j and tjj represent A and its uncertainty, respec- 
tively, for a given decay j, and s is a common scale pa- 
rameter introduced in the fit to account for a possible 
mis-estimate of &j. The convolution is defined by: 

S\(\j,(Tj) = J G(x - Xj,aj)exp ^—^j dx, (4) 

where As = ctb, and tb is the lifetime of the Af, (B°). 
The distribution of the uncertainty of A is modeled by an 
exponential function convoluted by a Gaussian. 

For the background, the mass distribution of the 
prompt component is assumed to follow a flat distribu- 
tion as observed in data when a cut of A >100 /im is 
applied. The non-prompt component is modeled with a 
second-order polynomial function. The A distribution is 
parametrized by the resolution function for the prompt 
component, and by the sum of negative and positive ex- 
ponential functions for the non-prompt component. A 
positive and a negative exponential functions model com- 
binatorial background, and an exponential function ac- 
counts for long-lived heavy flavor decays. The distribu- 
tion of the uncertainty of A is modeled by two exponential 
functions convoluted by a Gaussian. 

We minimize — 21n£ to extract: ct(A^) = 
365.1±g;* urn and ct(B°) = 450.0±^;^ /im. From 
the fits, we obtain ,s = 1.41 ± 0.05 for the Af, and 
s = 1.41 ± 0.03 for the B°. The numbers of signal decays 
are 171 ± 20 A b and 717 ± 38 B°. Figures [1] and [2] show 
the mass and A distributions for the Af, and B° candi- 
dates. Fit results are superimposed. Table |l] summa- 
rizes the systematic uncertainties in our measurements. 
The contribution from possible misalignment of the SMT 
detector was estimate to be 5.4 ^m 14]. We estimate 
the systematic uncertainty due to the models for the A 
and mass distributions by varying the paramctrizations 
of the different components: (i) the resolution function 
is modeled by two Gaussian functions instead of one, (ii) 
the exponential functions in the non-prompt background 
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FIG. 1: Invariant mass distribution for At (left) and B° (right) candidates, with the fit results superimposed. The inserts show 
the mass distributions after requiring X/a > 5. 




'■■ u-' u-J -U.l U U.I U.Z U.3 U.' 

Proper decay length (cm) Proper decay length (cm) 



FIG. 2: Proper decay length distribution for At (left) and B° (right) candidates, with the fit results superimposed. The shaded 
region represents the signal. 



are replaced by exponentials convoluted with the reso- 
lution function, (iii) a uniform background is added to 
account for outlier events (this has only a negligible ef- 
fect), (iv) the positive and negative exponentials describ- 
ing combinatorial non-prompt background are assumed 
to be symmetric, and (v) for the mass distribution of the 
non-prompt background, a linear function is used instead 
of the nominal quadratic form. To take into account cor- 
relations between the effects of the different models, a fit 
that combines all different model changes is performed. 
We quote the difference between the result of this fit and 
the nominal fit as the systematic uncertainty. 

The lifetime of the background events under the 
Ab(B°) signal is mostly modeled by events in the low 
and high mass sideband regions with respect to the peak. 
To estimate the effect of any difference between the life- 
time distributions of these two regions, we perform sep- 
arate fits to the A;, (B°) mass regions of 5.1-5.8 and 



5.4-6.1 GeV/c 2 (4.9-5.45 and 5.1-5.7 GeV/c 2 ) where the 
contributions from high and low mass background events 
are reduced, respectively. The largest difference between 
these fits and the nominal fit is quoted as the systematic 
uncertainty due to this source. 

We also study the contamination of the A& sample by 
B° events that pass the A& selection. From Monte Carlo 
studies, we estimate that 6.5% of B° events pass the A& 
selection criteria. However, the invariant mass of B° 
events which contaminate the Af, sample is distributed 
almost uniformly across the entire A& mass range, and 
their proper decay lengths therefore tend to be incorpo- 
rated in the long-lived component of the background. To 
estimate the effect due to this contamination, we remove 
from the A& sample any event which also passes the B° 
selection criteria, and we perform a fit to the remaining 
events. The difference between this and the nominal fit 
is quoted as the systematic uncertainty due to the con- 
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TABLE I: Summary of systematic uncertainties in the mea- 
surement of cr for At and B° and their ratio. The total uncer- 
tainties are determined by combining individual uncertainties 
in quadrature. 



Source 


A b (/mi) 


5° h 


Ratio 


Alignment 


5.4 


5.4 


0.002 


Distribution models 


6.6 


2.8 


0.020 


Long-lived components 


6.0 


13.6 


0.022 


Contamination 


7.2 




0.016 


Total 


12.7 


14.9 


0.034 



tamination. For the B°, we do not consider this source 
of systematic uncertainty since any event identified as A& 
is removed from the B° sample. 

We perform several cross-checks of the lifetime mea- 
surements. The J/tp vertex is used instead of the b 
hadron vertex, the mass windows are varied, the recon- 
structed b hadron mass is used instead of the PDG [3[ 
value, and the sample is split into different pseudora- 
pidity regions and different regions of azimuth. All re- 
sults obtained with these variations are consistent with 
our measurement. We also cross-check the fitting proce- 
dure and selection criteria by measuring the A& lifetime 
in Monte Carlo events. The lifetime obtained was con- 
sistent with the input value. 

The results of our measurement of the A;, and B° life- 
times are summarized as: 

cr(A fc ) = 365.1±ig;2, (stat) ± 12.7 (syst) /mi, (5) 
ct{B°) = 450.0^22:1 ( stat ) ± 14 - 9 (syst) /mi, 



from which we have: 



r(A b ) = 1.218^'^ (stat) ±0.042 (syst) ps, 
t{B°) = 1.501±g;g|l (stat) ±0.050 (syst) ps. 



(6) 



These can be combined to determine the ratio of life- 
times: 



r(i?°) 



0.811 



i-0.096 
-0.087 



(stat) ± 0.034 (syst), (7) 



where we determine the systematic uncertainty on the ra- 
tio by calculating the ratio for each systematic source and 
quoting the deviation in the ratio as the systematic un- 
certainty due to that source. We combine all systematics 
in quadrature as shown in Table [T] The main contribu- 
tion to the systematic uncertainty of the lifetime ratio is 
due to the long-lived component of the B° sample. This 
is expected since the B° is more likely than the A b to 
be contaminated by mis-reconstructed B mesons due to 
its lower mass. The ratio of lifetimes, using the world 



average B° lifetime t(B° 
r(A b ) 



1.527 ±0.008 ps [3j, is 



= 0.797 



+0.089 
-0.080- 



t(B°) ~ -0-080' ( 8 ) 

In conclusion, we have measured the Aj, lifetime in the 
fully reconstructed exclusive decay channel J/tp A. The 
measurement is consistent with the world average Q , and 
the ratio of A;, to B° lifetimes is consistent with the most 
recent theoretical predictions Q . 
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